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We propose a quantization procedure for the nucleon-scalar meson system, in which an arbitrary 
mean scalar meson field (j) is introduced. The equivalence of this procedure with the usual ones is 
proven for any given value of <j). By use of this procedure, the scalar meson field in the Walecka's 
RMFT and in Chin's RHA are quantized around the mean field. Its corrections on these theories 
are considered by perturbation up to the second order. The arbitrariness of (j) makes us free to fix 
it at any stage in the calculation. When we fix it in the way of Walecka's RMFT, the quantum 
corrections are big, and the result does not converge. When we fix it in the way of Chin's RHA, 
the quantum correction is negligibly small, and the convergence is excellent. It shows that RHA 
covers the leading part of quantum field theory for nuclear systems and is an excellent zeroth order 
approximation for further quantum corrections, while the Walecka's RMFT does not. We suggest to 
fix the parameter <j) at the end of the whole calculation by minimizing the total energy per-nucleon 
for the nuclear matter or the total energy for the finite nucleus, to make the quantized relativistic 
mean field theory (QRMFT) a variational method. 
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I. INTRODUCTION 

Relativistic Mean Field Theory (RMFT) is a fruit- 
ful and widely used theory in nuclear physics [1-14]. Its 
agreement with observational data is impressive. It seems 
to show that the nuclear data under certain energy scale, 
say several hundred MeV, may be roughly understood by 
hadron field theory. Of course, it should be a quantum 
theory of fields. But in RMFT, meson field operators are 
replaced by their expectation values, which are consid- 
ered as classical. RMFT is therefore semi-classical. A 
meaningful RMFT must be followed by quantum correc- 
tions due to meson field quantization, and the correction 
should not qualitatively change its agreement with obser- 
vational data. In this case, one may realize a quantum 
hadron field theory for nuclear systems under certain en- 
ergy scale, with RMFT to be the zeroth order contribu- 
tion, and calculate quantum corrections by perturbation. 

Quantum hadron field theory has already been de- 
veloped for nuclear physics in the usual loop expansion 
formalism [1-6], in which the mean meson field is the con- 
tribution of the tadpole diagrams. The attachment of 
nucleon loops (tadpole heads) on a nucleon line by the 
scalar meson lines (tadpole tails) changes the nucleon 
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mass. Since the tadpole head itself is also formed by a 
nucleon line, the RMFT calculation for a nuclear sys- 
tem means an infinite series of attachments of tadpoles 
on tadpoles. In this sense, RMFT is non-perturbative, 
like the method of self-consistent field widely used in 
quantum mechanics. A reasonable way of considering 
its quantum corrections is to quantize the field around 
its classical solution, like the consideration of residual 
interactions in many-body problems on the basis of its 
self-consistent field solutions, or like the consideration of 
quantum corrections in laser-atom interaction on the ba- 
sis of its semi-classical solutions [15]. This is a quantiza- 
tion procedure for nucleon field expanded on a complete 
set of single nucleon solutions in the mean meson fields 
(instead of in vacuum) and for meson fields around their 
suitably chosen non-zero mean values (instead of around 
the vacuum). This is a generalization of the usual quan- 
tization procedure around the vacuum, and is equivalent 
to it. Instead the loop expansion, we take the expan- 
sion scheme in which terms are classified according to 
the number of mesons in the intermediate states. This 
is something like the Tamm-Dancoff method[16, 17] used 
in particle physics. Numerical calculations here will be 
limited to the approximation in which only one meson 
appears in the intermediate state. This is in the spirit 
of the one boson exchange (OBE) idea in the traditional 
nuclear physics. This procedure relates with the RMFT 
and with the usual nuclear theory more close. Here we 
would check if and in which case the quantum correction 
on RMFT may be regarded as a small perturbation. 

Section II is a formulation of the quantization proce- 
dure for the scalar meson field around its mean value 
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and that for the nucleon field in mean scalar meson field. 
Its equivalence with the usual quantization procedure in 
vacuum is shown. In section III we apply this proce- 
dure to the (j-uj model[2] for nuclear matter, both the 
RMFT contribution and its quantum corrections are de- 
rived. Numerical results are given. Section V is devoted 
to conclusions and discussions. 



II. FORMULATION 

Consider a system consisting of a neutral scalar meson 
field $ and a nucleon field , its Lagrangian density in 
the nature units 7i = c = 1 is 

C = § (17^9^ - m) * + is^^a'^* - im^ $2 



with 



(1) 



Let us introduce an arbitrary classical mean value for 
the scalar meson field , which is assumed to be indepen- 
dent of space-time. Defining 

$' = $-.^, (2) 

we write the Lagrangian density (1) in the form 

C = * (i7^a^ - m') * + - ^m^ $'2 



with 



m = m-ga(, 
C ='^<^ 

9a 



(3) 

(4) 
(5) 



The quantization of the nucleon field in (1) and (3) 
seems to be the same. But the sets of eigenfunctions 
on which one expands the field operator \1/ are different. 
In the former case we expand ^' in terms of the com- 
plete set of eigenfunctions // (k) exp (ik ■ x) and 

—r^ — i>sT-(k) exp (— ik • x) of the single nucleon energy 

V(27r)3 

operator a - (— iV) + /3m in vacuum: 



*(x) = j -^5][c«,(k)«,.(k)e"'-^ 



*(x) 



+rft^(k)^;,,(k)e-i''-^ 
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^ [ds^(k)i;k 



+4^(k)?xksre" 



-ik-x"| 



(6) 



(7) 



s and T are spin and isospin indices respectively, the nu- 
cleon spinor states Ust (k) and Vst (k) satisfy equations 

(a-k-h^m)u,^(k) =a;(fc)Mg^(k), (8) 
(a • k - 13m) vsr (k) = uj {k)vsr (k) , (9) 



Lj{k) = + m2 , (10) 

ul^{k)Us'r'{k) = Sss'Srr' (H) 
vl{k)v,,r'0^)=S,s'Srr', (12) 
ul{k)v,,,,ik) = vl{k)u,,r'{k) = , (13) 

quantization conditions are 

Cs7-(k)Cs/^/(k') + Cs'T'Oi')Csr(k.) 
= 4r(k)(is'T'(k') + ds'r'(k')dsr{'k) 
= Csr(k)4'r'(k') + 4'r'(k')CsT(k) 

= c«,(k)4,,(k') +4,,(k')c,.(k) =0, (14) 
c«,(k)cl,,,(k') + ct,,,(k')c,,(k) 

= d,,(k)4,^,(k') +4,,(k')rf«r(k) 

= 6ss'Srr'S{k--k'). (15) 

The vacuum state |0) is defined to be the eigenstate of 
annihilation operators with zero eigenvalues. It means 



c,,(k)|0) =d,,(k)|0) =0. 



(16) 



To make the; expcictation value bo zero in vacuum, the 
Hamiltonian density is expressed in terms of normal 
products, in which annihilation operators Csr(k) and 
(isr(k) always stand on the right of creation operators 
Cg,^,(k') and (i^,^,(k'). The nucleon sector of the Hamil- 
tonian density is therefore 

■Hi = : 'H'l {-id ■ V + (3m) - ga'^^'i : 
= {-id ■ V -I- (3m) * - Qa^^'^ 



m 



\/FT 



(17) 



products sandwiched between two colons are defined to 
be normal. For nuclear matter, ^ — A (neutrons and 
protons with spin up and down), and for neutron matter, 
7 = 2. 

In the later case, it is in the classical scalar meson field 
0, we may quantize \E' in the same way, but have to substi- 
tute m',u)', c^^(k), (ig^(k) and |0)' for to, w, Csx(k), (isr(k) 
and |0) respectively. The nucleon sector of the Hamilto- 
nian density is therefore 

H[ = :*^(-ia-V-F/3m')*-5<7*'**:' 
= {-id ■ V + /3m) * - fifcr^^* 
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TO 



k^dk. 



(18) 



Products sandwiched between : and :' are defined to be 

normal in the sense that annihilation operators c^^(k) 
and rfs^(k) are on the right of creation operators Cy^,(k') 
and<V,(k'). 
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H 



The difference between these two expressions gives 

(x/PT^ - Vfc^ + k'^dk 



Jo 

Jo 



m 



m 



l2 



m 



/2 



(19) 



The integrals have been worked out analytically. Deleting 
a fourth degree polynomial in (j) and by the renormal- 
ization of scalar meson field energy density, we obtain 



A = 



167r2 



i4 1 

mm- 



(20) 

7 

m^{m — m') — -m^{m — m')^ 



13 , . . Q 

+ —m{m — m ) 



25 
12 



/\4 



(m — m') 



(21) 
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+ m (m — m') — -m(m — m') 



+ — m-m'f 

D 



The usual quantization of the scalar meson field $ and 
its canonically conjugated variable 11 in vacuum is per- 
formed by the expansions 



4>(x) 



d^k 



n(x) 



V(27r)32a;<,(fc) 
+at(k)e-'''-"] , 

d^fc 



[a(k), 



,ik-x 



a(k)e*-^ 



^i^i [at(k)e-"'-^ 



with 



(23) 

(24) 

(25) 

(26) 
(27) 

Instead, one may quantize the scalar meson field $' 
around a classical value (p by the same procedure. It is to 
substitute 11' and a'(k) for 11 and a(k) in equations 
(23), (24), (26), and (27). The relation (2) gives 



Quantization conditions are 

a(k)a(k') - a(kXk) =0, 
a{k)a\k') - at(k')a(k) = S{k - k') 



a(k) = a'(k) + , 



(28) 



(29) 



Since the difference of $ and is an additive c-number, 
we have 



(V$)^ n'^ + (V$') 



S'\2 ./ 



(31) 



in which the normal products sandwiched in : : are de- 
fined in terms of annihilation operators a(k) and creation 
operators a^(k), while the normal products sandwiched 
in : :' are defined in terms of annihilation operators a'(k) 
and creation operators a'^{\i). The Hamiltonian of the 
nucleon-scalar meson system is therefore 

H = j {■■'^^ (-i<5 • V + I3m) * - g^^^^ : 
= j {-.'^^-ia-V + I3m')^ :' -A 



+g„ [{B + C) : ^^'-^ :'] } d^x , 



(32) 

This equation shows the equivalence of the quantization 
procedure in and around the mean meson field with the 
usual one in and around the vacuum. Write 



H = + H' , 



(33) 



(22) with 



= /a;'(fc)5^[c;Uk)c;,(k)+<t(k)<,(k)]d3fc 

ST 

Lj„{k)a"'{k)a'{k)d^k 



1 



-m 



2 J,2 



- ^ d^a; 



uj'{k) = Vfc^ + m'- 



= 9.j [( 



B + C)<^'- : :'] d^x. 



(34) 

(35) 
(36) 



The unperturbed Hamiltonian i?*^"^ depends only on the 
mean value (j) of the scalar meson field, no single scalar 
meson appear in it. It is therefore the RMFT Hamilto- 
nian for the nucleon-scalar meson system. Since A is the 
change of the nucleon vacuum energy due to the appear- 
ance of (/), the term containing it represents a quantum 
effect. The usual RMFT discards this term, but Chin's 
RHA takes it. RHA is therefore an extended RMFT, or 
shortly the ERMFT. The perturbation H' contains the 
quantum correction due to the fluctuation of the scalar 
meson field around its mean value </> and the nucleon- 
scalar meson interaction. The theory considering both 
and H' is therefore a quantized RMFT, or shortly 
the QRMFT. 



III. QUANTUM CORRECTIONS OF THE 
WALECKA a-LO MODEL FOR NUCLEAR 
MATTER 

The zeroth order ground state of an iso-symmetric 
static uniform nuclear matter is defined by the eigen- 

state \rj) of TJ*-"-* with nucleons filled in the positive en- 
ergy Fermi sea from the bottom up to the Fermi surface 



of momentum 77, and without any single scalar meson. 
The contribution of H' is considered by perturbation up 
to the second order. 

A vector meson field is needed to stabilize the nu- 
clear matter. The Walecka cr-w model treats the nuclear 
matter as a nucleon-scalax meson-vector meson field sys- 
tem by RMFT. a is an iso-singlet scalar meson with a 
mass about several hundred MeV to be determined by 
the comparison of the theoretical calculation with nu- 
clear data, while w is an iso-singlet vector meson with a 
mass of 783MeV. The contribution of the vector meson 
on the energy density of a static uniform nuclear mat- 
ter is gi^rP' /2m^, in which is the w-nucleon coupling 
constant, muj is the uj meson mass, and n = 777^/677^ is 
the nucleon number density of nuclear matter. This con- 
tribution comes from a second order perturbation of the 
nucleon-vector meson coupling. The energy pcr-nuclcon 
of the nuclear matter in units ofc = ?i = m= lis then 



e = ei -I- £2 -I- 



9Z: ^ 

2m},, n 



A gl{n-B-CY 



in which 



X_ 
2772 




(37) 



(38) 



is the average energy of a nucleon in the nuclear matter, 
with X = m'/m = 1 — g4>/m to be the efi'ective nucleon 
mass in unit of its free mass m; 



£2 



W9a 



(39) 



is the energy of the mean scalar meson field per nucleon; 



n = 



7X 
47r2 



l\/rfT)(^ - In 




(40) 



a--o 



FIG. 1: Feynman diagram 



TABLE I: Paramotors and calculatcxl properties of iso- 
symmetric nuclear matt(>r at e(iuilil)riuiii density 





X 


A' (MeV) 


RMFT 362.6 278.0 


0.538 


554 


RHA 230.3 149.0 


0.730 


456 


QRMFT 230.3 149.0 


0.730 


456 



due to the appearance of the mean field 0, and therefore 
shows the renormalized vacuum fluctuation. The sum of 
the first four terms has been considered by Chin in his 
relativistic Hartree approximation (RHA) [7]. As we ex- 
plained before, it is an extended version of the RMFT. 
The last term is our new, it comes from the second order 
perturbation of H', and is somewhat complex. There 
is a square of a sum with three terms. Expanding the 
square, one obtains nine terms. The one with comes 
from the interaction between nucleons in positive energy 
Fermi sea by the exchange of a scalar meson in medium, 
it is the quantum of the field This is an eff'ect of 
the in-medium OBEP, and is shown in FIG.l. The term 
with comes from the change of interaction energy be- 
tween negative energy nucleons in vacuum by exchang- 
ing a scalar meson in medium, due to the appearance 
of the mean field </>. The term with comes from the 
change of the scalar meson field vacuum. Other terms 
show the mixing and interference of these three effects. 
Altogether, they show typical quantum effects. The sum 
of these five terms is the result of QRMFT in the second 
order perturbation approximation. 



is the scalar nucleon number density of the nuclear mat- 
ter; with 



A = 



B = 



C = 



7 



167r2 
13 



X^lnx+(1-X)-^(1-X)' 



+y(l-x)- 



25 4 

Y^(i-x) 



7 
4772 

11 

"6 

ml{l - x) 



X'lnx + (1-X)-^(1-X)' 



+^(i-x)^ 



(41) 

(42) 
(43) 



The sum of the first three terms on the right of (37) is 
exactly the usual RMFT result considered by Walecka[2] , 

while the last two terms arc quantum corrections. The 
fourth term comes from the change of vacuum energy 



IV. NUMERICAL RESULTS 

The energy e per-nucleon is a function of x and 77. x 
is determined by the condition 



dx 



= 



(44) 



of the energy minimization. It makes e be a function 
of 77 alone, and therefore be a function €{n) of nucleon 
number density n. This is the nuclear equation of state. 
The model parameters are reduced to two independent 
dimensionless parameters = m'^gl/m'^ and = 
m?g^/m^. They are chosen to reproduce the binding 
energy b = 15.986MeV per-nucleon and the equilibrium 
density 3/ (47rrQ) of nuclear matter with rg — 1.175fm[18]. 

The resulting parameters are listed in the TABLE I. 
The first line is given by Walecka's RMFT, which is ob- 
tained by approximating e by the sum of its first three 
terms in (37) before to be substituted into the variational 



1.5 

n//i 



0.6 

«/«„ 



FIG. 2: Equation of state for iso-symmetric uniform nuclear 
matter, given by RMFT (line 1) and by RHA or QRMFT 
(line 2) 

condition (44). The second line is given by Chin's RHA, 
which is obtained by approximating e by the sum of its 
first four terms in (37) before to be substituted into the 
variational condition (44). The third line is given by our 
QRMFT, which is obtained by substituting the whole 
expression (37) for e into the variational condition (44). 
The sets of parameters are adjusted to reproduce the 
same set of nuclear data listed above. We see that quan- 
tum corrections notably change the parameters. The last 
two columns show two of the most important data calcu- 
lated by corresponding methods. The low density equa- 
tion of state is almost determined by the compression 
modulus K[8]. One of the shortcomings of the Walecka's 
RMFT is that it gives a too large K value. Various meth- 
ods, such as the introduction of non-linear terms [8], have 
been introduced to overcome this shortcoming. How- 
ever we see that it is noticeably reduced by the quantum 
corrections, without introducing any additional parame- 
ter. Another important property is the effective mass x, 
which determines the energy dependence of the optical 
potential. The value 0.730 in last two lines is consistent 
with the reasonable range 0.65 < x < 0.75[19], while that 
of the Walecka's RMFT is too small. 

Looking at TABLE I, one amazedly sees that the re- 
sults of RHA and QRMFT are the same. This is also 
shown in FIG. 2. The equations of state for RHA and 
QRMFT are the same. It means that the quantum cor- 
rection on RHA given by the last term of (37) is practi- 
cally zero, and therefore the RHA solution of the problem 
defined by the Lagrangian (1) is an excellent zeroth order 
approximation for the next order quantum correction. 

On the other hand. Fig. 3 shows that the quantum 
corrections on RMFT given by the last two terms of (37) 
are quite big, and make the equation of state qualitatively 
distorted. The series of quantum corrections seems not 



FIG. 3: Line 1: The RMFT equation of state. Line 2: The 
equation of state obtained by adding the fourth term of (37) 
on line 1. Line 3: The equation of state obtained by adding 
the fifth term of (37) on line 2. 

converge. It means that the Walecka's RMFT solution of 
the problem defined by the Lagrangian (1) is not a good 
zeroth order approximation for the next order quantum 
corrections. 

Of course, these results are for the problem defined 
by the Lagrangian (1). The situation should be checked 
case by case. The smallness of the correction on RHA is 
due to the cancelation between n, B and C. In the prob- 
lem considered above, the cancelation is almost complete. 
However, the completeness will not be always. There is 
evidence showing that if the Lagrangian (1) is general- 
ized in a way with the term ^m^^'^ being substituted 
by 

C/($) = im2$2 (i + ai$ + a2$2) , (45) 

the quantum corrections on RHA may be non-zero but 
still small, for non-zero ai and 02. 

V. CONCLUSIONS AND DISCUSSIONS 

Shortly speaking, our conclusions are: 

1. The RHA covers the leading part of quantum field 
theory for nuclear systems, but the simple RMFT 
does not. 

2. The residual interaction H' in Eq. (36) may be con- 
sidered by perturbation theory, with the expansion 
scheme used in traditional nuclear physics. 

3. QRMFT, which we proposed here, is therefore an 
easy and practical way for handling nuclear prob- 
lems until quark degrees of freedom become impor- 
tant. 
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In our formalism, the mean field (j) is a free parameter. 
The equivalence (32) is true for its any value. This makes 
us have right to choose its value so that the solution of 
our problem becomes easy and accurate. RMFT, RHA 
and QRMFT chose it in different ways, and then obtained 
different accuracies and convergence properties for their 
solutions. Among them the QRMFT is the best. In 
this way we solve the problem of non-convergence, which 
the loop expansion procedure also encountered before [6]. 
Different value of tj) means different vacuum in the nu- 
clear matter. The theory with a change of vacuum must 
be non-perturbative. The vacuum |0)' in nuclear matter 
satisfies the condition 

a'(fc)|0)' = 0. (46) 

By (28) we see 

a(fc)|O)' = 0/^3M!53(k)|o)'. (47) 



It means |0)' is a coherent state of the original scalar 
mesons. Scalar mesons develop a Bose condensation in 
nuclear matter. 



Beside the scalar meson, any kinds of mesons may de- 
velop Bose condensation in nuclear matter under suitable 
condition. The most often talked is the 7r-condensation. 
We have even proven that the 7r-condcnsation may de- 
velop in the Walecka model[20, 21]. In non-uniform nu- 
clear systems, for example the finite nuclei, the vector 
meson condensation may develop. In these cases, the 
quantum fluctuation of the condensed mesons should also 
be considered. The meson condensation is described by 
a classical meson field, and the quantum fluctuation is 
described by its quantum corrections. The method de- 
veloped here may be useful for considering these correc- 
tions. 
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